Universal conductivity and dimensional crossover in multi-layer graphene 
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We show, by exact Renormalization Group methods, that in multi-layer graphene the dimensional 
crossover energy scale is decreased by the intra-layer interaction, and that for temperatures and 
frequencies greater than such scale the conductivity is close to the one of a stack of independent 
layers up to small corrections. 
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Recent experiments [H have found that the optical con- 
ductivity in multi-layer graphene is essentially constant 
in a wide range of frequencies and equal to Nao, with 

2 

N < 5 the number of layers and do = remarkably, 
a value both independent from the intra-layer hopping 
t ^ 3eV and inter-layer hopping ^ O.SeV and de- 
pending only on the conductivity quantum h/e'^ and the 
number of layers N. In the absence of interaction, when 
= 1 (monolayer graphene) this observation finds a nice 
explanation in the effective description of fermions on the 
honeycomb lattice in terms of massless Dirac fermions in 
2 + 1 dimensions, whose conductivity in the limit of zero 
frequency is exactly ctoj as was shown in [ij. The band 
structure of multi-layer graphene is completely different 
and the description in terms of Dirac fermions breaks 
down at low energies; in the case of bilayer graphene, 
for instance, 0], a quadratic spectrum instead of a lin- 
ear relativistic one is found for certain kind of hopping 
terms. However as pointed out in P, 0] in multi- layer 
graphene there is a dimensional crossover scale t± and at 
energy scales greater than t± the conductivity of multi- 
layer graphene is equal to the one of a stack of N inde- 
pendent layers up to small corrections. 

The above reasoning neglects the role of the interac- 
tions among charge carriers, which are well known to be 
rather strong and producing observed important effects 
in several physical observables [H']. In the case of mono- 
layer graphene there are several papers discussing the role 
of interactions on the conductivity , but in the case 

of multi-layer this problem is much less studied, even if 
interactions in multi-layer graphene are believed to rad- 
ically alter the low energy properties and spontaneous 
symmetry breaking is expected 



ll|,ll2|. 



In the present paper we will consider a system of sev- 
eral layers in which the electrons can hop from a plane 
to another with coupling t-^, and t"*" << t; each layer is 
described in terms of interacting fermions on the honey- 
comb lattice with an electromagnetic (e.m.) interaction 
isl 14 1 . We will show that the interaction produces a 
renormalization of the dimensional crossover scale and 
the conductivity of multi-layer graphene is equal to the 
one of a stack of iV independent layers up to small cor- 



rections at energy scales greater than 



^ t ' 



(1) 



where 77 > is the exponent of the wave function renor- 
malization. Note that the crossover scale is decreased 
by the intra-layer interaction, a phenomenon strongly re- 
sembling what happens in fermionic chains, see e.g. [l5| . 
Even if our results are found for values of the parameters 
close to the infrared fixed point of monolayer graphene 
(that is, the Fermi velocity large enough), they provide 
an evidence that the many-body interaction preserves 
and even enforces the conductivity universal behavior, 
in qualitative agreement with experiments flj. 

We consider the Hamiltonian for multi-layer graphene 



N 



(2) 



a=l 



where Ha describe a single graphene layer; more exactly 
fermions hopping on the honeycomb lattice interacting 
through an e.m. field introduced via the Peierls substi- 
tution 
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Ha^-t V at b- ^ e'^^^o ^"■^•^=(^+''^3.o)'i« -^c.c. 



xGAa 



(3) 



dk ipi e 

kA,a 



for 



electrons with plane index a = 1,2 and sitting at the 
sites of the two triangular sublattices and of a 
honeycomb lattice. We assume that A^ has basis vectors 
fi 2 = i(3, ±%/3) and that Ab = Aa+Sj, with Si = (1, 0) 
and (52,3 — ^{"^,^VS) the nearest neighbor vectors; B 
is the first Brillouin zone and A^. and (/)q, are the respec- 
tively the vector field and the coulomb potential on the 
plane a. 

Finally H± describes the hopping between graphene 
layers and is assumed of the form 



N-l 



(4) 
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where ^ is a bilinear operator describing the hopping 
from layer a to layer a + 1, see e.g. The planar 

paramagnetic current is defined as 



N 

Jp=iet^ 



where rj^ 



.3 — 1-e' 



.pi,- 



ipSj 



'InUai b- -. -6+ . ) 

J lp\ x,a 3j^Sj,a x+5,,a a;, a/ 



The two components of the para- 



magnetic current Jp will be seen as the spatial compo- 
nents of a "space-time" three-components vector Jp-,^, 
/i = 0, 1, 2, with Jpfl = epp and pp the density operator. 

If = e^f'^^'Os^e-^o^A, with x = (xo,f), we 
denote by (OxV ■ ' ' Ox""*) the thermodynamic limit of 

S-iTr{e-''^T(04V ' ' ' where S = Tr{e"''^} and 

T is the operator of fermionic time ordering. The 
current-current functions Kf^^{Tp) is defined as the 2D 
Fourier transforms of (Jx./^; Jy,v) p and the conductivity 
is Q (here l,m = 1,2): 



--4=-[i<'i.™(c^,0)-ii:i,™(0,0)] (5) 
3v3 w 



where 3-\/3/2 is the area of the hexagonal cell of the 
honeycomb lattice. In our notations, p = (aj,p), with 
a; € the Matsubara frequency. 

It is convenient to introduce the following Grassman 
integral 



(6) 



where: tp^ are Grassman variables (k — (fco,fc)) and 
P{dip) is the fermionic gaussian integration with inverse 
propagator 



vD,{k) ikf) 



(7) 



2^ ^uy., ~ ^T,j=12 3^'''^'^'^^^'' (note that 



with w = ft rj(fc) 

(7(k) is singular only at the Fermi points k 
(0, ^,±-^)); Hp, = 0, 1,2, A^{p) are gaussian variables 

with propagator (5a./3i5pi/^^^ where x ^'Ct an an ultravio- 
let cut-off; finally V is the interaction whose explicit form 
can be easily inherited from H. By suitable derivatives 
with respect to J the current-current correlation can be 
obtained; note that we have exploited gauge invariance 
to write the photon propagator in the Feynman gauge. 
The generating function ^ can be computed by exact 
Rcnormalization Group methods. After the integration 
of the fields A^^\ A^''+^'> we get 

(8) 



where P{dijA~^^) is the fermionic integration with prop- 
agator, if r = ± is the valley index 

Xh{^) ( iko Vhn*{k'+p^)Y 

where is effective wave rcnormalization and Vh 
the effective Fermi velocity, P{dA'^~'^^) is the gauge 
field integration with propagator Sa^/^S^^i, ^2 \ p\^ ' with 
X;i(k'), x;i(p) smooth cut-off functions with support 
smaller than t2'^; moreover V is the effective potential 
which has the form 



(9) 



where the kernels W^m depends on the effective charge 
at scale h e^^h and the effective hopping t^. We have now 
to describe the integration of the field ^l)^^\A^'^^ and in 
this way we will iteratively define the effective constants 
Zh: VhjCf^^h, tj[. In order to do that we have to decompose 

^(h) y(h) ^ ^ ^yih) ^^^^ U = 1 - C] CV'^''^ IS 

the relevant or marginal part of the effective interaction 
while TZV^'^^ is the irrelevant part. Generally this decom- 
position is dictated by the naive scaling dimension which 
in the present case is given by 



D 



n — m 



(10) 



C should select the terms with positive or vanishing di- 
mension D. However, if /i^ is the temperature scale 
fit — 2^'*'^ and if the temperature verifies the condition 



t2'^^ > ti^ 



(11) 



where is the hopping at scale /i, there is an improve- 
ment with respect to naive power counting, and certain 
terms which are dimensionally relevant or marginal are 
indeed irrelevant. In order to verify this fact, we can split 
the kernels as 



(12) 



where W^m^^ is obtained from Wn% setting t-^ ~ 0. 
Note that W^m^^ in correspondence of external fields 
with different plane index are vanishing. 
We define the C operator as 



£W'£^(k') = V[^2\"i^"^(0) 
so that 

Ah)/-,!-. _ rTl/(i)('»),'u'\ _ w('^)('')/ 



r{a){h). 



(13) 



W2^7(k') = [W^27^'^^(k') - W!fl'-'''m + M^fi)"')(k') 

(14) 

The first term in the r.h.s. of (ITil) can be rewritten as 
k' ■dW^)}i^\ and this produces an improvement ~ 2'' ~^ 
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in the bound of the kernel, if h' is the scale of the momen- 
tum, which is sufficient to make it irrelevant. Similarly 
the second term in (fM)) . namely W2^l^^\'k'), has an ex- 
tra t^2-'' < 2('''3-'*)(i-''), r] = 0(e'2) (see below) with 

respect to the bound for ^^2^]^ , which again is enough to 
make it irrelevant; therefore, the true marginal contribu- 
tion is given by the r.h.s. of (fT3)) . Regarding the terms 
quadratic in the gauge fields, 

>C<2^ (P) = WI^T^ (0) + pdwl.f''^ (0) (15) 



where we have used that wlf^'"'^^ (0) has an extra {2~^tj^ 
with respect to the naive dimension. Finally the terms 
quadratic in the fermionic variables, if they have the same 
plane index then 



£W'2^_';)(k') = w^27^"\o) -I- k'aw27^"^(k' 



,(a){h) 



{a){h). 



(16) 



where we have used that in W^l'^^^^ there is an extra gain 
0((t^2~'')^, due to the conservation of the plane index 
a. On the other hand for the quadratic terms with ex- 
ternal fields corresponding to j-^ we define £W^2 Q^(k') = 
M^2g\o), while for the other terms with different layer 
index£W'2^o^(k') = 0. 

The terms in d^wi";^^^'^ and diW!^'^^'^^'^ with both fields 
with the same plane index are included in the free 
fermionic intergration and produces the new effective 
wave function renormalization Zh~i and Fermi velocity 
Vfi-i', therefore we get after rescaling 



„v('^'(v'z;~rv''-''',A(^'')-Hj)) 



(17) 



with 



a=l 



where J^^p'' is the intra-layer current and j^a'^'^ the 
inter-layer current. Note that by construction the effec- 
tive Fermi velocity Vh, effective wave function renormal- 
ization and effective charge e-^ = e^ ^ Si^e the same 

by ex 

we get, see [iJl Vh = O(e^) and 



as m the = case and by exploiting Ward Identities 



Zh 



= l + 0(e) 



Z 



ZhVh 



= l + 0(e) 



(19) 



Moreover, see [14|, the wave function renormalization di- 
verges with a power law with a critical exponent and the 



effective Fermi velocity increases up to the light velocity 
with a power law 



with r/ 



+ 0{el^) and 



2e' 



hh 



(20) 



+ 0(e^-oo)- 



Regarding the flow of tf^ we obtain 



Zh 



(21) 



with \I3[^\ < Cii^e*^(t^2-'*)2. It is easy to see by induc- 
tion that \Zhtj^ — i"*"! < C2t'^e^. Assume indeed that it 
is true for k > h: therefore 



\ttiZH-i~t^\<2t^Ce'^J2^ti2 



-k\2 



(22) 



k=h 



from which the inductive assumption follows. Note that 
the effective hopping, even if relevant in the RG sense 
according to naive power counting, remains small in this 
region of temperatures. Moreover, from (jlip we obtain 
the condition between the temperature and the hopping 



>t^(^)T^ = r 



(23) 



As the flow of the effective parameters corresponding to 
the relevant and marginal operators is bounded, the fol- 
lowing bound is obtained, for h > hjj (order by order in 
the renormalized expansion) 



^ j dxK'),(x)| < C2''(3- 



(24) 



We apply the above bound to the conductivity, which is 
given by 



2 1 



3%/3w 



da;o(e'' 



l)i^,,,(x) (25) 



We can decompose ^^^^(x) = K\^{'x.) + K\1{-k) where 

K'fl{^) is obtained from Ki^i{-x) setting — Q and, for 
any M 



\K, 



(a) 



(x)| < r^'^ 



Zh' l + (2''|x|)A^ 



(2"|x|) 



M 



(26) 



The above estimates can be derived from the dimensional 
bound (|24|) : roughly speaking there is, with respect to 
([24l) . an extra 2'^'' due to a lacking integration and a 

which can be extracted from 



decaying factor i+(2fe|x|)M 
the chain of propagators connecting the external fields. 
In the bound for ^^^'^''(x) there are also two extra 
Therefore 



(27) 



4 



where 



2 1 



3%/3 w 
so that for t*-^ < 



\R^\u;)\ < C (i^)" J 

h—hs 



mu^ (28) 



2(2+2r;)/i 



(2''|x|)5 



Therefore for i* -'- < (3^^ « lo the conductivity is the 
one of a stack of independent layers plus a negligible cor- 
rections, in qualitative agreement with the observations 
in[l|. 
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